Molecular control of cellulosic fin morphogenesis in ascidians.

BACKGROUND: The tunicates form a group of filter-feeding marine animals closely related to vertebrates. They share with them a number of features such as a notochord and a dorsal neural tube in the tadpole larvae of ascidians, one of the three groups that make tunicates. However, a number of typical chordate characters have been lost in different branches of tunicates, a diverse and fast-evolving phylum. Consequently, the tunic, a sort of exoskeleton made of extracellular material including cellulose secreted by the epidermis, is the unifying character defining the tunicate phylum. In the larva of ascidians, the tunic differentiates in the tail into a median fin (with dorsal and ventral extended blades) and a caudal fin. RESULTS: Here we have performed experiments in the ascidian Phallusia mammillata to address the molecular control of tunic 3D morphogenesis. We have demonstrated that the tail epidermis medio-lateral patterning essential for peripheral nervous system specification also controls tunic elongation into fins. More specifically, when tail epidermis midline identity was abolished by BMP signaling inhibition, or CRISPR/Cas9 inactivation of the transcription factor coding genes Msx or Klf1/2/4/17, median fin did not form. We postulated that this genetic program should regulate effectors of tunic secretion. We thus analyzed the expression and regulation in different ascidian species of two genes acquired by horizontal gene transfer (HGT) from bacteria, CesA coding for a cellulose synthase and Gh6 coding for a cellulase. We have uncovered an unexpected dynamic history of these genes in tunicates and high levels of variability in gene expression and regulation among ascidians. Although, in Phallusia, Gh6 has a regionalized expression in the epidermis compatible with an involvement in fin elongation, our functional studies indicate a minor function during caudal fin formation only. CONCLUSIONS: Our study constitutes an important step in the study of the integration of HGT-acquired genes into developmental networks and a cellulose-based morphogenesis of extracellular material in animals.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis.

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Effects of color variation and physiological state on ascidian microbiomes.

Ascidians, known for their color variation, host species-specific microbial symbiont communities. Some ascidians can also transition into a nonfiltering (resting) physiological state. Recent studies suggest that the microbial symbiont communities may vary across different physiological states and color morphs of the host. The colonial ascidian, Polyclinum constellatum, which exhibits several color morphs in the Caribbean Sea, periodically ceases its filtering activity. To investigate if color variation in P. constellatum is indicative of sibling speciation, we sequenced fragments of the ribosomal 18S rRNA and the mitochondrial cytochrome oxidase subunit I genes. Additionally, we sequenced a fragment of the 16S rRNA gene to characterize the microbial communities of two common color morphs (red and green) in colonies that were either actively filtering (active) or nonfiltering (resting). Phylogenetic analyses of both ascidian genes resulted in well-supported monophyletic clades encompassing all color variants of P. constellatum. Interestingly, no significant differences were observed among the microbial communities of the green and red morphs, suggesting that color variation in this species is a result of intraspecific variation. However, the host's physiological state significantly influenced the microbial community structure. Nonfiltering (resting) colonies hosted higher relative abundances of Kiloniella (Alphaproteobacteria) and Fangia (Gammaproteobacteria), while filtering colonies hosted more Reugeria (Alphaproteobacteria) and Endozoicomonas (Gammaproteobacteria). This study demonstrates that microbial symbiont communities serve as reliable indicators of the taxonomic state of their host and are strongly influenced by the host's feeding condition.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids.

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi.

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Specification of distinct cell types in a sensory-adhesive organ important for metamorphosis in tunicate larvae.

The papillae of tunicate larvae contribute sensory, adhesive, and metamorphosis-regulating functions that are crucial for the biphasic lifestyle of these marine, non-vertebrate chordates. We have identified additional molecular markers for at least 5 distinct cell types in the papillae of the model tunicate Ciona, allowing us to further study the development of these organs. Using tissue-specific CRISPR/Cas9-mediated mutagenesis and other molecular perturbations, we reveal the roles of key transcription factors and signaling pathways that are important for patterning the papilla territory into a highly organized array of different cell types and shapes. We further test the contributions of different transcription factors and cell types to the production of the adhesive glue that allows for larval attachment during settlement, and to the processes of tail retraction and body rotation during metamorphosis. With this study, we continue working towards connecting gene regulation to cellular functions that control the developmental transition between the motile larva and sessile adult of Ciona.

Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs.

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

Genome-wide analysis of early vascular tunic repair and regeneration for Botrylloides digenesis reveals striking similarities to human wound healing.

The early stages of regeneration after injury are similar to those of wound healing. The ascidian Botrylloides diegensis can regenerate an entire adult from a small fragment of vascular tunic following the removal of all zooids in an injury-induced regeneration model. We investigated the molecular and cellular changes following injury to determine the differences between the healing process and the initiation of whole-body regeneration (WBR). We conducted transcriptome analysis at specific time points during regeneration and wound healing to identify differentially expressed genes (DEGs) and the unique biological processes associated with each state. Our findings revealed 296 DEGs at 10 h post-injury (hpi), with 71 highly expressed in healed tissue and 225 expressed during the WBR process. These DEGs were predicted to play roles in tissue reorganization, integrin signaling, extracellular matrix organization, and the innate immune system. Pathway analysis of the upregulated genes in the healed tunic indicated functional enrichment related to tissue repair, as has been observed in other species. Additionally, we examined the cell types in the tunic and ampullae in both tissue states using histology and in situ hybridization for six genes identified by transcriptome analysis. We observed strong mRNA expression in cells within the WBR tunic, and in small RNA-positive granules near the tunic edge. We hypothesized that many of these genes function in the compaction of the ampullae tunic, which is a pivotal process for WBR and dormancy in B. diegensis, and in an immune response. These findings establish surprising similarities between ascidian regeneration and human wound healing, emphasizing the potential for future investigations into human regenerative and repair mechanisms. This study provides valuable insights into the gene sets specifically activated during regeneration compared to wound healing, shedding light on the divergent activities of these processes.

Drop it all: extraction-free detection of targeted marine species through optimized direct droplet digital PCR.

Molecular biomonitoring programs increasingly use environmental DNA (eDNA) for detecting targeted species such as marine non-indigenous species (NIS) or endangered species. However, the current molecular detection workflow is cumbersome and time-demanding, and thereby can hinder management efforts and restrict the "opportunity window" for rapid management responses. Here, we describe a direct droplet digital PCR (direct-ddPCR) approach to detect species-specific free-floating extra-cellular eDNA (free-eDNA) signals, i.e., detection of species-specific eDNA without the need for filtration or DNA extraction, with seawater samples. This first proof-of-concept aquarium study was conducted with three distinct marine species: the Mediterranean fanworm Sabella spallanzanii, the ascidian clubbed tunicate Styela clava, and the brown bryozoan Bugula neritina to evaluate the detectability of free-eDNA in seawater. The detectability of targeted free-eDNA was assessed by directly analysing aquarium marine water samples using an optimized species-specific ddPCR assay. The results demonstrated the consistent detection of S. spallanzanii and B. neritina free-eDNA when these organisms were present in high abundance. Once organisms were removed, the free-eDNA signal exponentially declined, noting that free-eDNA persisted between 24-72 h. Results indicate that organism biomass, specimen characteristics (e.g., stress and viability), and species-specific biological differences may influence free-eDNA detectability. This study represents the first step in assessing the feasibility of direct-ddPCR technology for the detection of marine species. Our results provide information that could aid in the development of new technology, such as a field development of ddPCR systems, which could allow for automated continuous monitoring of targeted marine species, enabling point-of-need detection and rapid management responses.

Salpa genome and developmental transcriptome analyses reveal molecular flexibility enabling reproductive success in a rapidly changing environment.

Ocean warming favors pelagic tunicates, such as salps, that exhibit increasingly frequent and rapid population blooms, impacting trophic dynamics and composition and human marine-dependent activities. Salp blooms are a result of their successful reproductive life history, alternating seasonally between asexual and sexual protogynous (i.e. sequential) hermaphroditic stages. While predicting future salp bloom frequency and intensity relies on an understanding of the transitions during the sexual stage from female through parturition and subsequent sex change to male, these transitions have not been explored at the molecular level. Here we report the development of the first complete genome of S. thompsoni and the North Atlantic sister species S. aspera. Genome and comparative analyses reveal an abundance of repeats and G-quadruplex (G4) motifs, a highly stable secondary structure, distributed throughout both salp genomes, a feature shared with other tunicates that perform alternating sexual-asexual reproductive strategies. Transcriptional analyses across sexual reproductive stages for S. thompsoni revealed genes associated with male sex differentiation and spermatogenesis are expressed as early as birth and before parturition, inconsistent with previous descriptions of sequential sexual differentiation in salps. Our findings suggest salp are poised for reproductive success at birth, increasing the potential for bloom formation as ocean temperatures rise.

Stem cell-mediated development, regeneration, chimerism, and aging in the colonial chordate Botryllus schlosseri.

Stem cells are units of biological organization, responsible for tissue and organ development and regeneration. I study stem cell biology, aging, and the evolution of immunity using the colonial chordate Botryllus schlosseri as a model system. This organism is uniquely suited for this study because it is closely related to vertebrates, undergoes weekly cycles of stem cell mediated regeneration, is long lived and has a recognition system and robust immune system. I have led the Botryllus genome project and developed a novel method to obtain a synthetic long read sequence, identified Botryllus stem cells and stem cell niches, isolated the gene that controls self/non self-recognition and characterized its immune system on the cellular and molecular levels. Recently, I led the Botryllus atlas project to characterize the two developmental pathways, embryogenesis (sexual) and blastogenesis (asexual), revealing the unique molecular landscapes for each developmental mode and investigated the molecular clock and neurodegeneration pathways in young and old colonies and investigated the molecular clock and neurodegeneration pathways in young and old colonies. These results and the resources we developed are used by my lab and others to further study stem cell and immune cell properties during development, regeneration, transplantation, and aging.

De novo genome assembly and comparative genomics for the colonial ascidian Botrylloides violaceus.

Ascidians have the potential to reveal fundamental biological insights related to coloniality, regeneration, immune function, and the evolution of these traits. This study implements a hybrid assembly technique to produce a genome assembly and annotation for the botryllid ascidian, Botrylloides violaceus. A hybrid genome assembly was produced using Illumina, Inc. short and Oxford Nanopore Technologies long-read sequencing technologies. The resulting assembly is comprised of 831 contigs, has a total length of 121 Mbp, N50 of 1 Mbp, and a BUSCO score of 96.1%. Genome annotation identified 13 K protein-coding genes. Comparative genomic analysis with other tunicates reveals patterns of conservation and divergence within orthologous gene families even among closely related species. Characterization of the Wnt gene family, encoding signaling ligands involved in development and regeneration, reveals conserved patterns of subfamily presence and gene copy number among botryllids. This supports the use of genomic data from nonmodel organisms in the investigation of biological phenomena.